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Abstract 

The Stern diatomic sequence is closely linked to continued fractions via the Gauss 
map on the unit interval, which in turn can be understood via systematic subdivi¬ 
sions of the unit interval. Higher dimensional analogues of continued fractions, called 
multidimensional continued fractions, can be produced through various subdivisions 
of a triangle. We define triangle partition-Stern sequences (TRIP-Stern sequences for 
short) from certain triangle divisions developed earlier by the authors. These sequences 
are higher-dimensional generalizations of the Stern diatomic sequence. We then prove 
several combinatorial results about TRIP-Stern sequences, many of which give rise to 
well-known sequences. We finish by generalizing TRIP-Stern sequences and presenting 
analogous results for these generalizations. 


1 Introduction 

Stern’s diatomic sequence (defined in Section 2.2) stems from the study of continued fractions 
and has a number of remarkable combinatorial properties, as seen in Northshield [45]. There 
are many different multidimensional continued fraction algorithms, and they serve a num¬ 
ber of different purposes ranging from simultaneous Diophantine approximation problems 
(see Lagarias [37]) to attempts to understand algebraic numbers via periodicity conditions 
(see the third author’s [22]) to automata theory (see Fogg [19]). This paper concerns a 
generalization of Stern’s diatomic sequence dehned using triangle partition maps, a family 
of multidimensional continued fractions that includes most of the well-known multidimen¬ 
sional continued fractions presented in Schweiger [49]. For background on multidimensional 
continued fractions, see Schweiger [49] and Karpenkov [29]. 

For background on the properties of Stern’s diatomic sequence, see Lehmer [38]. For back¬ 
ground on Stern’s diatomic sequence in the context of continued fractions, see Northshield 
[45]. Knauf found connections between Stern’s diatomic sequence and statistical mechanics 
[32, 33, 34, 35, 36], though Knauf called the sequences Pascal with memory, which is a more 
apt description. The connection between Stern’s sequence and statistical mechanics was fur¬ 
ther developed in Contucci and Knauf [8], Esposti, Isola and Knauf [15], Fiala and Kleban 
[17], Fiala, Kleban and Ozluk [18], Garrity [20], Guerra and Knauf [24], Kallies, Ozluk, Peter 
and Syder[28], Kleban and Ozluk [31], Mayer [41], Mendes France and Tenenbaum [42, 43], 
Prellberg, Fiala and Kleban [47], and Prellberg and Slawny [48]. Other important earlier 
work was done by Allouche and Shallit [2, 3], who showed that Stern’s sequence is 2-regular. 

There seems to have been little work on extending Stern’s diatomic sequence to mul¬ 
tidimensional continued fraction algorithms. The hrst generalization of Stern’s diatomic 
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sequence was for a type of multidimensional continued fraction called the Farey map, in 
Garrity [21], The Farey map is not one of the multidimensional continued fractions that we 
will be considering. Another generalization used the Monkemeyer map, in Goldberg [23]. 

This paper uses the family of multidimensional continued fractions called triangle par¬ 
tition maps (TRIP maps for short) [10] to construct analogous sequences. As mentioned, 
many, if not most, known multidimensional continued fraction algorithms can be put into 
the language of triangle partition maps; thus, the collection of TRIP maps is a rich family. 
In Section 2, we give a quick overview of continued fractions. Stern’s diatomic sequence and 
how the two are related. Section 3 reviews triangle partition maps and triangle partition 
sequences. Section 4 introduces the construction of TRIP-Stern sequences. In Section 5, we 
give a more pictorial description of TRIP-Stern sequences. Section 6 contains results about 
the maximum terms and locations thereof for each level of the TRIP-Stern tree. Section 7 
discusses minimum terms and locations thereof. Section 8 examines sums of levels of the 
sequence. Section 9 determines which lattice points appear in the TRIP-Stern sequence for 
the triangle map, a multidimensional fraction algorithm discussed below. Section 10 intro¬ 
duces a generalization of the original TRIP-Stern sequence. Finally, we close in Section 11 
with some of the many questions that remain. 


2 Continued fractions and Stern’s diatomic sequence 

Nothing in this section is new. In the hrst subsection, we review continued fractions in 
order to motivate, in part, the dehnition of triangle partition maps given in Section 3. In 
the second subsection, we review the classical Stern’s diatomic sequence and show how it is 
linked to continued fractions. This link is what this paper generalizes. 


2.1 Continued fractions and subdivisions of the unit interval 


All of the content in this subsection is well-known. 

Let a be a real number in the unit interval I = (0,1]. The Gauss map is the function 
G : (0,1] —)■ [0,1) dehned by 


G{a) = - 


a 


where [xj denotes the floor function, meaning the greatest integer less than or equal to x. 
Subdivide the unit interval into subintervals 


h. = 


k + G 

for k a positive integer. If a 6 /*,, then the Gauss map is simply G{a) = xhe continued 
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fraction expansion of a is 


a = 


Oq 


ai- 


“2 + - 


where a E lag, G{a) G G{G{a)) G Ia2, ■ ■ ■■ (If «, under the iterations of G, is ever zero, 
then the algorithm stops.) 

We now want to translate the dehnition of the Gauss map into the language of two-by-two 
matrices, which can be more easily generalized. Set 


Vi = 


0 


and V 2 = 


1 y ^ ^1 

We have the standard identihcation of a vector in to a real number via 

X 


-)■ 


y 


provided of course that |/ 7 ^ 0. Then we think of the two-by-two matrix 

0 1 


V = (Vi,V2) = 


1 1 


as being identihed to the unit interval I. Set 

0 1 


Fn = 


1 1 


and Fi = 


1 1 
0 1 


Then, by a calculation, we have 


VF^-^Fo = 


1 1 
k k + 1 


which can be identihed to the subinterval Ik- Further, by a calculation, we have that 


V{VFt^Fo)-^{ " 


-)■ 


-k 1 
1 0 

—ka + 1 
a 

1 — ka 


a 


and thus captures the Gauss map. 


5 






Using the matrices Fq and Fi, we can also interpret the Gauss map as a method of 
systematically subdividing the unit interval. This interpretation leads to the classical Stern 
diatomic sequence. Note that 


VFq = {vi,V2)Fo 

= {vi,V2) ( J j ) • 


= {V2,VI+V2) 



and 


VFi = {vi,V2)Fi 

= J). 


= {vi,Vi+V2) 



We can interpret VFq as the half interval (1/2,1) and VFi as the half interval (0,1/2). If 
we iterate multiplying by Fq and Fi, then we get the following at the next step: 


Each real number a E I can be described by a sequence (foGiU 2 , • • •) of zeros and ones, 
where, for all n, the number a lies in the subinterval coming from VFi^Fi^Fi^ ■ ■ ■ F^^. (We 
are being somewhat sloppy with issues of a being on the boundaries of these subintervals. 
Such issues do not affect what is going on.) We can link the sequence (foUiU 2 , • • •) with 
a’s continued fraction expansion as follows. Let 1^ denote a sequence of k ones. Then our 
sequence (fo, H, * 2 , • • •) can be written as 

(*0UiU2,...) = (l'^°,0,l"b0,l"U0,...), 


with each kj a non-negative integer. (It is important that we allow a kj to be zero.) Then 
we have 

1 


a = 


/cn T 1 T 


fcl+l-l- 


w 


^2 + 1 + — 
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For example, the sequence (1,1, 0, 0,1,1,1, 0,...) can be written as (1^, 0,1°, 0,1^, 0,...), 
and we have 

1 

a =-^-. 

2 + 1 +^ 


4 +- 


Thus, continued fractions can be interpreted as a systematic method for subdividing an 
interval using two-by-two matrices. Multi-dimensional continued fractions, as we will see, 
are systematic subdivisions of triangles determined by three-by-three matrices. 


2.2 Stern’s diatomic sequence 

In this section, we will briefly review Stern’s diatomic sequence (number A002487 in Sloane’s 
Online Encyclopedia of Integer Sequences). In particular, we highlight the link between 
Stern’s diatomic sequence and continued fractions. The classical Stern’s diatomic sequence 
oi, 02 , Os,... is the sequence dehned by oi = 1 and, for n > 1, 


CL2n 

^2n+l “1“ ^n+1* 


Stern’s diatomic sequence is linked to the Stern-Brocot array, which is an array of frac¬ 
tions in lowest terms that contains all rationals in the interval [0,1]. Starting with the 
fractions ^ and j on the level, we construct the level by rewriting the {n — 1)®* level 

with the addition of the mediant between consecutive pairs of fractions from the {n — 1)®* 
level. Here, the mediant of two fractions x and § refers to the fraction In the Stern- 
Brocot array, the mediant of consecutive fractions is always in lowest terms. Below are levels 
0 through 3 of the Stern-Brocot array: 


0 


1 ’ 







0 




1 



1 ’ 




2 ’ 



0 


1 


1 


2 

1 ’ 


3 ’ 


2 ’ 


3 ’ 

0 

1 

1 

2 

1 

3 

2 

1 ’ 

4 ’ 

3 ’ 

5 ’ 

2 ’ 

5 ’ 

3 ’ 


1 

1 

1 

1 

1 

1 

1 

1 


The denominators of the Stern-Brocot array form Stern’s diatomic sequence. Many of 
the combinatorial properties of this sequence are presented in Northshield [45]. 

The hrst row of the array can be thought of as either the unit interval, or, as above, the 

two-by-two matrix V = ( ^ j V The second row can be thought of as the two subintervals. 


(0,1/2) and (1/2,1), or, as the two matrices VFq and VFi. Similarly, the third row gives 
us four subintervals, each corresponding to one of the matrices VFqFq^VFqFi^VFiFq and 
VFiFi. The pattern continues. 
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To be more precise, let Sn,k denote the fraction in the level of the Stern-Brocot 
tree. One can use the Stern-Brocot array to express the continued fraction expansion of a 
real number in [0,1] as follows: let a G [0,1]. The 1®* level of the Stern-Brocot tree divides 
the unit interval in two as the sub intervals [0, |) and [|, 0]. Label the hrst interval 0 and the 
second interval 1. The 2'^'^ level divides the unit interval into four sub intervals: [0, |), 

[|, |), and [|, y]- Label these intervals 00, 01, 10, and 11 respectively. The level divides 
the unit interval into [0, s„^i),..., !]• We label the interval [sn,k, Sn^k+i) with a sequence 

of O’s and I’s, where the hrst n — 1 digits mark the label of the interval containing [sn^k, Sn,k+i) 
on the (n — 1)®* level, and where the last digit is 0 or 1 depending on whether [sn^k, Sn,k+i) is 
in the left or right half of that interval, respectively. Recording the inhnite sequence of O’s 
and Ts that corresponds to any number a in [0,1] yields a sequence encoding the continued 
fraction expansion of a, as in described in Northshield [45]. Thus, Stern’s sequence is linked 
to subdivisions of the unit interval. Our generalizations of Stern’s sequence will be linked to 
subdivisions of a triangle. 


3 Review of triangle partition maps 

TRIP-Stern sequences can be interpreted geometrically in terms of subdivisions of a triangle. 
(This section closely follows Sections 2 and 3 from Dasaratha et ah [10].) In this section, we 
describe the triangle division and triangle function, as dehned in Garrity [22] and further 
developed in Chen et ah [7] and Messaoudi et al. [44]. We then discuss how “permutations” 
of this triangle division generate a family of multidimensional continued fractions called 
triangle partition maps (TRIP maps for short) - which were introduced in Dasaratha et al. 
[10, 11] - and studied in Jensen [27] and in Amburg [4]. This will give us the needed notation 
to dehne TRIP-Stern sequences in the next section. 


3.1 The triangle division 


The triangle division generalizes the method for computing continued fractions via the above 
systematic subdivision of the unit interval from the previous section. Instead of dividing the 
unit interval, we now use a 2-simplex, i.e., a triangle. As discussed in earlier papers, this 
triangle division is just one of many generalizations of the continued fraction algorithm. 
Dehne 

A* = {{bo, bi, 62) ■ bo > bi > b 2 > 0}. 


The set A* can be thought of as a “triangle” in using the projection map vr : ^ 

dehned by 


7r(6o,&i,&2) 


(b_, bA 
Wbo)' 


The image of A* under tt. 


A = {{l,x,y) : 1 >x >y > 0}, 


is a triangle in with vertices (0, 0), (1, 0), and (1,1). Thus vr maps the vectors 


= 0 1 ^2 = 1 , ^^3 = 1 


0 


0 


to the vertices of A. The change of basis matrix from triangle coordinates to the standard 
basis is 

'1 1 A 


(ui ^2 W3) = I 0 1 1 I . 

.0 0 1 , 


Now consider the matrices 


0 0 1 


1 0 1 


Ao = I 1 0 0 and /li = 0 1 0 

0 11 / \ 0 0 1 

Applying Aq and Ai to (wi V2 v^) yields 

(wi V2 V3)Ao = {v2 V3 vi + V3) and (wi V2 ^3)^1 = (wi V2 vi + V3). 


This gives a disjoint bipartition of A under the map vr, as seen in the diagram below. 

( 1 . 1 ) 



This is the hrst step of the triangle division algorithm. Now consider the result of applying 
Ai to the original vertices of the triangle k times followed by applying Aq once. The vertices 
of the original triangle A are thus mapped as follows: 


Afc = {(1, x, y) e A : 1 — a: — /c|/ > 0 > 1 — X — (/c + 1)|/}. 

We let T : A^ —)■ A denote a collection of maps, where each is the following bijection 
between the subtriangle A^ and A: For any point (x,?/) G A^ under the standard basis, 
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first change the basis to triangle coordinates by multiplying the coordinates by (ui V 2 V 3 ), 
apply the inverse of Aq and the inverse of A\, and then hnally change the basis back to the 
standard basis. That is, 

T{x,y) = 7r((l X y) ((ui V 2 n 3 )(^o)"^(^i)"''(wi V 2 ^ 3 )"^)^). (1) 

Rewriting yields the following: 

\X X J 

where k = I I. 

L jy J 

This triangle map is analogous to the Gauss map for continued fractions. Using T, dehne 
the triangle sequence, {tk)k>o, for an element {a, /5) G A by setting tn to be the non-negative 
integer satisfying G Thus, 

{a, (3) G AtQ,T(a, (3) G A^^, T(T(a, (3)) G ... 


3.2 Incorporating permutations 

The previously described triangle division partitions the triangle with vertices Vi,V 2 , and ^3 
into triangles with vertices V 2 ,v^, and Vi + and Vi,V 2 , and Vi -|- ^ 3 . This process assigns 
a particular ordering of vertices to the vertices of the original triangle and the vertices of 
the two triangles produced. By considering all possible permutations we generate a family 
of 216 maps, each corresponding to a partition of A. 

Specihcally, we allow a permutation of the vertices of the initial triangle as well as a 
permutation of the vertices of the triangles obtained after applying A^ and Ai. First, we 
permute the vertices by a G S '3 before applying either Aq or Ai. Once we apply either Aq or 
Ai, we then permute by either tq G S 3 or ti G 5*3, respectively. This leads to the following 
dehnition: 

Definition 1. For every (a, ro,ri) G Ff, dehne 

Fq = aAoTo and Fi = aAiTi 


by thinking of a, Tq, and Ti as column permutation matrices. 


In particular, we dehne the permutation matrices as follows: e = 


0 1 0 

100 ), (13) = 
0 0 1 


0 0 1 \ / 1 0 0 

010 , (23) = 0 0 1 I , (123) 

10 0 / \ 0 1 0 


(132) = 


1 0 0 \ 
010 , ( 12 ) = 
0 0 1 / 

(0 1 0 \ 

( 0 0 1 1 , and 

Vl 0 0 ; 


0 0 1 
1 0 0 
0 1 0 


10 



Note that applying Fq and Fi partitions any triangle into two subtriangles. Thus, for 
any (a, ro,ri) G S'!, we can partition A using the matrices Fq and Fi instead of Aq and Ai. 
This produces a map that is similar to, but not identical to, the triangle map from section 
3.1. We call each of these maps a triangle partition map, or TRIP map for short. Because 
IS'sl^ = 216, the family of triangle partition maps has 216 elements. 

One of the main goals of Dasaratha et ah [10] is showing that this class of triangle partition 
maps includes well-studied algorithms such as the Monkemeyer map, and in combination, 
these triangle partition maps can be used to produce many other known algorithms, such as 
the Brun, Parry-Daniels, Guting, and fully subtractive algorithms. 

3.3 TRIP sequences 

We use the triangle division for a particular set of permutations to produce the corresponding 
triangle sequence. Recall the “subtriangle A^.” from the original triangle map. We generalize 
the dehnition of A^. as follows: for any (a, Tq, Ti), let A^ be the image of the triangle A under 
F^Fq. We now define functions T : A —)■ A mapping each subtriangle A^ bijectively to A. 
First, let us formalize our subtriangles. 

Definition 2. Let Fq and Fi be generated from some triplet of permutations. We define 
An to be the triangle with vertices given by the columns of (ui V 2 v^)FA'Fq. 

We are now ready to define triangle partition maps. 

Definition 3. We define the triangle partition map T'(T,ro,ri by 

Ta,ro,ri(a:,2/) = vr((l,x,|/)((ni V2 V3,)Fq^FR^{ vi V2 

when (x, y) G A^. 

We take the transpose of the matrix (ui V 2 V'i)FQ^FR^{yi V 2 because our matrices 

have vertices as columns but they are multiplied by the row vector {l,x,y). 

This definition facilitates the following: 

Definition 4. For any an (a,/3) G A, define to be the non-negative integer such that 
[To-,ro,Ti]”(Q:,/^) is in A*^. The TRIP seguence of {a, (I) with respect to (a, ro,ri) is {tk)k>o- 

4 TRIP-Stern sequences 

4.1 Construction of TRIP-Stern sequences 

Definition 5. For any permutations (a, ro,ri) in S 3 , the triangle partition-Stern seguence 
{TRIP-Stern seguence for short) of (a, ro,ri) is the unique sequence such that ai = (1,1,1) 
and, for n > 1, 
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®2n Ojn ■ -^0) 
®2n+l ~ ■ Fl- 


The level of the TRIP-Stern sequence is the set of with 2” ^ < m < 2"'. 

Each choice of {a, Tq, Ti) produces some TRIP-Stern sequence. The level terms of the 
TRIP-Stern sequence give the hrst coordinate of the vertices of the subtriangles of A after 
n divisions. These terms are the denominators of the convergents of the triangle partition 
map defined by (a, tq, ti). Thus, TRIP-Stern sequences can be used to test when a sequence 
of triangle subdivisions converges to a unique point. This involves a simple definition and 
restatement of Theorem 7.2 in Dasaratha et al. [10]. 

Definition 6. For each n-tuple v = (ii,..., A) of O’s and I’s, dehne 




For a TRIP-Stern term a„, the length of the binary representation v such that A(n) = an 
gives the level of an- If there are multiple binary representations of an, then choose the 
representation that gives the term of the sequence. 

4.2 Examples 

We start by examining the TRIP-Stern sequence for the triangle map, which is given by the 
identity permutations (a, Tq, Ti) = (e, e, e). The first few terms of the sequence {an)n>i are 



We arrange these as follows: 


ai = (1,1,1) 



02 = ( 1 , 1 , 2 ) 


03 = ( 1 , 1 , 2 ) 




or as 


04 —(1,2,3) 05 —(1,1,3) Og —(1,2,3) 07 —(1,1,3) 

A 



A( 0 ) 


A(l) 




A( 00 ) 


A( 01 ) A( 10 ) A(ll) 

12 


Note that because the two triples on level 1 are the same, the left and right subtrees are 
symmetric. To see the connection to the triangle division, recall the matrix 


(■^1 V2 Us) 



whose columns are the three vertices of A. Applying Aq and Ai to (ui V 2 V 3 ) yields 

/I 1 2\ 

and (ui V 2 ^ 3)^1 = I 0 1 1 . 

Vo 0 1 ) 

Repeating this process yields matrices at each step n whose columns give the vertices of 
the subtriangles at the rA^ division of A. The figure below shows this subdivision. Addition¬ 
ally, the top row of each of these step matrices gives a term in the n^^Aievel TRIP-Stern 
sequence for (e, e,e). 


(I 1 2 \ 

(Ul V2 ^3)^0 =1 1 1 

Vo 1 1/ 



For another example, consider the permutations (a, ro,ri) = (13,132,132). Recall here 
that we divide the triangle using the matrices Fq = 13Aol32 and Fi = 13Ail32. 

Applying Fq and Fi to (ui V 2 vA) yields 


and 


/I 1 0\ /I 2 1 

(ui V2 vA)Fo = (ui ^2 t^3) I 0 0 1 j = I 0 1 1 

Vo 1 0/ Vo 1 0 

/O 1 0\ /I 2 1\ 

(ui V2 vAjFi = (ui ^2 '^ 3 ) 1 0 0 = 1 1 1 

Vo 1 1/ Vo 1 V 
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These permutations give the well-studied Monkemeyer map, as described in Dasaratha 
et ah [10]. Here, the first few terms are 

( 1 , 1 , 1 ), ( 1 , 2 , 1 ), ( 1 , 2 , 1 ), ( 1 , 2 , 2 ), ( 2 , 2 , 1 ), ( 1 , 2 , 2 ), ( 2 , 2 , 1 ),... 

Similarly, these terms now give the hrst coordinate of the vertices of the subtriangles of 
A according to the Monkemeyer algorithm. The following figure shows this division. 



We will examine TRIP-Stern sequences generated by each of the 216 maps. As we will 
see in Lemma 9, many properties of each of these TRIP-Stern sequences can be captured by 
only examining the 36 sequences associated with permutation triplets of the form (e, ro,ri) 
for ro,ri G S 3 . The following table, which illustrates the behavior of Fq and Fi for maps 
of the form (e, ro,ri) on the initial seed {a,b,c) (which equals (1,1,1) in the case of regular 
TRIP-Stern sequences), will come in handy throughout the paper. 
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{e,To,Ti) 

(e,e,e) 

(e,e, 12 ) 

(e,e, 13) 

(e,e,23) 

(e,e,123) 

(e,e,132) 

(e, 12 , e) 
(e, 12 , 12 ) 
(e,12,13) 
(e,12,23) 

(e, 12,123) 
(e, 12,132) 
(e, 13, e) 
(e,13,12) 
(e,13,13) 
(e,13,23) 

(e, 13,123) 

(e, 13,132) 

(e,23,e) 

(e,23,12) 

(e,23,13) 

(e,23,23) 

(e, 23,123) 
(e,23,132) 
(e,123,e) 

(e, 123,12) 
(e, 123,13) 
(e, 123,23) 
(e, 123,123) 
(e, 123,132) 
(e,132,e) 

(e, 132,12) 
(e, 132,13) 
(e, 132,23) 
(e, 132,123) 
(e, 132,132) 


{a,b,c)Fo 
( 6 , c, a + c) 
(fe, c, a + c) 
( 6 , c, a + c) 
(fe, c, a + c) 
( 6 , c, a + c) 
(fe, c, a + c) 
(c, 6 , a + c) 
(c, 6 , a + c) 
{c,b,a + c) 
(c, 6 , a + c) 
{c,b,a + c) 
{c,b,a + c) 
(a + c, c, b) 
(a + c, c, b) 
(a + c, c, b) 
(a + c, c, b) 
(a + c, c, b) 
(a + c, c, b) 
( 6 , a + c, c) 
( 6 , a + c, c) 
(fe, a + c, c) 
( 6 , a + c, c) 
(fe, a + c, c) 
( 6 , a + c, c) 
(a + c, 6 , c) 
(a + c, 6 , c) 
(a + c, 6 , c) 
(a + c, 6 , c) 
(a + c, 6, c) 
(a + c, 6 , c) 
(c, a + c, 6 ) 
(c, a + c, 6 ) 
(c, a + c, 6 ) 
(c, a + c,b) 
(c, a + c, 6 ) 
(c, a + c, 6 ) 


(g, b, c)Fi 
(g, 6 , g + c) 
( 6 , g, g + c) 
(g + c, b, a) 
(g, g + c, 6 ) 
(g + c, g, 6 ) 
( 6 , g + c, g) 
(g, b,a + c) 
{b, g, g + c) 
(g + c, b, a) 
(g, g + c, 6 ) 
(g + c, g, 6 ) 
( 6 , g + c, g) 
(g, b,a + c) 
{b, g, g + c) 
(g + c, b, a) 
(g, g + c, 6 ) 
(g + c, g, 6 ) 
{b, g + c, g) 
(g, 6 , g + c) 
( 6 , g, g + c) 
(g + c, b, a) 
(g, g + c, 6 ) 
(g + c, g, 6 ) 
( 6 , g + c, g) 
(g, b,a + c) 
{b, g, g + c) 
(g + c, b, a) 
(g, g + c, 6 ) 
(g + c, g, b) 
( 6 , g + c, g) 
(g, 6 , g + c) 
( 6 , g, g + c) 
(g + c, b, a) 
(g, g + c, 6 ) 
(g + c, g, 6 ) 
( 6 , g + c, g) 


Table 1; Behavior of Fq and Fi 
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4.3 Alternative definitions 

An alternative way to define the TRIP-Stern sequences is in terms of matrix generating 
functions. The advantage of this method is that it allows us to connect terms in a TRIP- 
Stern sequence with the product of Fi matrices that produced the given sequence and to 
dehne the sequence non-recursively. 

Let 

P{x) = Fo + Fix, 

where x commutes with Fq and Fi. Now it is clear that any integer can be uniquely expressed 
as a sum of the form 2 "" + k, where 0 < /c < 2 "’. 

Definition 7. We dehne the TRIP-Stern sequence for (a, tq, ti) to be the unique sequence 
dehned by 

O-2'^+k = ( 1 , 1 , 1 ) • B, 

where B is the coefhcient of x^ in the product P{x)P{x^) ■ ■ ■P(x"'). Then a 2 "+A: is the 
term of the level of the TRIP-Stern sequence. 

For instance, take (a, ro,ri) = (e,e, e). Then 

P{x')P{x^') = (Aq -I- Aix)(Ao -I- Aix^') = Aq AiAqx AqAix'^ A'^x^ 

and so the terms of the 2 "'^ level are given by 

ilAA)Al (l,l,l)AiAo, (l,l,l)AoAi, and ( 1 , 1 ,1)4?. 

5 A more pictorial approach 

We need the technical dehnitions from the last section for the proofs given in the rest of the 
paper and in order to use Mathematica to discover many of our formulas. But there is a 
more intuitive approach to generate particular TRIP-Stern sequences via subdivisions of a 
triangle. 

Let us hrst examine the classical Stern diatomic sequence analog. Start with an interval 
/ whose endpoints are a and b. For Stern’s diatomic sequence, we set a = 6 = 1. As before, 
we set / = (a, b). 

I I 


a b 1 1 


We then subdivide the interval and add together the weights of the endpoints, getting two 
new intervals, which by an abuse of notation we write as J(0) = (a, a-\-b) and /(I) = {a-\-b, b). 
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/(i) 

- 1 

a a + b b 


h 


+ 


1 2 


H 
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We can continue, getting four new subintervals: 


/(OO) /(ll) 

1 2 a + b a + 2b i 
I-1-1-1-1 

a \a + b \ b 

/(Ol) /(lO) 


/(OO) 


t 2 

/(Ol) /(lO) 


/(ll) 

p t , 

t 1 


The classical Stern diatomic sequence corresponds to I, /(O), /(I), -f(OO), /(Ol), /(lO), /(ll),... 

We now see how to get the analogous geometric picture for TRIP maps. We will concen¬ 
trate on the TRIP-Stern sequence for the triple (e, e, e). Similar pictures, though, will work 
for any {a,To,Ti) G 

For (e, e, e) we have Fq = Aq, Fi = Ai. For any vector (a, b, c), we know that 

(a, b, c)Aq = (6, c, a -|- c) and (a, b, c)Ai = (a, b,a + c). 

We will think of {a,b,c) as vertices of a triangle A. By a slight abuse of notation, we will 
write A = (a,6, c). In the diagrams for this section, the ones on the left side are for the 
general case of any a, b, and c, while the right side is in the special case when a = 6 = c = 1. 
We have 



We let A(0) be the triangle with vertices (a, b, c)Ao = (/, c, a-l-c) and A(l) be the triangle 
with vertices (a, b, c)Ai = (a, b^a + c). Pictorially, we have 



17 


















We let A(00) be the triangle with vertices (a, b, c)AqAq = {c, a + c, a + b + c), A(01) be 
the triangle with vertices (a, b, c)AqAi = (6, c, a + 6 + c), A(10) be the triangle with vertices 
{a,b,c)AiAQ = {b,a + c,2a + c) and A(ll) be the triangle with vertices {a,b,c)AiAi = 
(a, b, 2a + c) 

Pictorially, we have 


c 



1 



A(10) 

The TRIP-Stern sequence for 


A(10) 

e,e,e) is simply A, A(0), A(l), A(00), A(01), A(10), 


A(ll), ... 

We just add the appropriate vertices to generate new subdivisions. This is in direct 
analog to the classical Stern diatomic sequence. 

As another example, let us look at the triangle partition Stern sequence for (12, e, e). We 
hrst need the matrices Fq and Fii 


and 


/ 0 1 0 
Fo = (12)Ao =10 0 

\ 0 0 1 

/ 0 1 0 
Fi = (12)Ai =10 0 

\ 0 0 1 


Now for any vector (a, b, c), we have 


0 0 1 
1 0 0 
0 1 1 

1 0 1 
0 1 0 
0 0 1 


1 0 0 \ 
0 0 1 
0 11 / 

0 1 0 \ 
10 1 
0 0 1 / 


(a, b, c)Fo = (a, c,b + c) and (a, b, c)Fi = {b, a,b + c). 


We still let A have vertices a, b, and c, but now the subtriangle A(0) will have vertices a, c, 
and b + c and the subtriangle A(l) will have vertices a, b, and b + c. Pictorially, we have 
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Let us do one more iteration. We have that A(0, 0) will have vertices {a,b + c,b + 2c), 
A(0,1) will have vertices (c, a, b + 2c), A(l, 0) will have vertices {b, b + c, a + b + c) and A(l, 1) 
will have vertices {a,b,a + b + c) since 


(a, b, c)FoFo 


(a, b, c)FoFi 


(a, b, c)FiFo 


(a, c,b + c)Fo 
{a,b + c,b + 2c), 

(a, c,b + c)Fi 
(c, a,b + 2c), 

(6, a,b + c)Fo 

(6, 6 + c, a + 6 + c). 


and 

(a, b, c)FiFi = {b, a,b + c)Fi 

= {a,b,a + b + c). 

By continuing these subdivisions, we get the TRIP-Stern sequence for (12,e, e). Pictori- 
ally, we now have 


c 



1 



6 Maximum terms and positions thereof 


This section examines the maximum terms in every level of a TRIP-Stern sequence, as 
well as the positions of those maximum terms within the level. Recall for an n-tuple v = 
(R,..., in) of O’s and I’s that A(n) = (1,1, l)Fi-^Fi^ ■ ■ ■ Fi^, which can be written as A(n) = 
{bi{v),b 2 {v),b^{y)). Let |n| denote the number of entries in v. 

Definition 8. The maximum entry on level n of a TRIP-Stern sequence is 

= max max bjiv). 

\v\=n ie{l,2,3} ^ ^ 

Thus, for example, the sequence (m„)„>o for the TRIP-Stern sequence for (e, e, e) begins 

(1,2,3,4,6,9,13,19,28,...). 
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We can easily extend this kind of analysis to all 216 TRIP-Stern sequences. Numerically, 
it appears that there are only eight possible row maxima sequences for all 216 TRIP-Stern 
sequences. The following lemma will be needed before presenting our results on maximum 
row sequences. 

Lemma 9. Let (a, ro,ri) G Rf have row maxima sequence {mn)n>o and suppose n G S 3 . 
Then (kct, roK“^, also has row maxima sequence (m„)n>o- 

Proof. This follows by direct calculation, since for any v = (R,..., i„), 

= ( 1 , 1 , 1 ) • KaAi^Ti^K~^ ■ KaAi^Ti^K~^ ■ ■ ■ KaAi^Ti„K~^ 

= (1, 1, 1) ■ cr^iiPi • crA2p2 • • ■ 

A(cr,ro,ri) 

Since is just a permutation, it follows that the maximal component of A(Ko-,ToK-i,TiK-i)('y) 
is the same as that of A(o-,T-o,Ti)('y)- CH 

Recall from Section 4 that the {n+ 1)®* level of a TRIP-Stern sequence is generated from 
the by applying Fq or Fi to the sequence at the mA" term at level n. Denote by Um the 
position of this term. In other words, a TRIP-Stern sequence is a binary tree with exactly 2 
children at each node rim, with the hrst (left) child generated by applying Fq to the triplet 
at and the second generated by applying Fi. We will show that there exist three distinct 
paths through the trees containing all the maximal terms. We will then present recurrence 
relations for select row maxima sequences. 

First, we show that there exist three distinct paths through the trees generated by triangle 
partition maps of the form (e,ro,ri) containing the maximal terms. For each of the three 
classes of path, we will write out an explicit proof for one map generating such path; proofs 
for the rest of the maps are similar simple calculations and will be omitted. 

Theorem 10. There exists a path through each of the trees generated by 26 triangle partition 
maps of the form (e, ro,ri) that contains the maximal terms. Namely, these paths are as 
follows: 

1. For the eleven TRIP-Stern sequences (e, e, 13),(e, 13,13),(e, 13,123), (e, 23,13), (e, 23,123), 
(e, 23,132), (e, 123,13), (e, 123,123), (e, 132,13), (e, 132,123) and (e, 132,132), always 
select the right edge. 

2. For the twelve TRIP-Stern sequences (e, e, e), (e, e, 12), (e, e, 23), (e, e, 123), (e, e, 132), 

(e, 12,e), (e,12,12), (e,12,13), (e,12,23), (e, 12,123), (e, 12,132) and (e, 132, 23), al¬ 
ways select the left edge. 

3. For the three TRIP sequences (e, 13,12), (e, 23, 23) and (e, 123,e), alternate between 
first selecting the left edge and then the right edge. 
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Proof. 1. Always select the right edge. 

Consider (e, 13,123). Then {xi,X 2 ,X 3 )Fq = {xi + X 3 ,X 3 ,X 2 ) and {xi,X 2 ,X 3 )Fi = (xi + 
X 3 ,Xi,X 2 ). Since at the zeroth step we start with {xi,X 2 ,x^) = (1,1,1), it is clear that 
the rightmost triplet will contain the maximal term at each level n since this will lead 
to the greatest rates of growth for each of X 2 „ and x^^ as n increases. 

Similar arguments can be made for (e, e, 13), (e, 13,13), (e, 23,13), (e, 23,123), (e, 23,132), 
(e, 123,13), (e, 123,123), (e, 132,13), (e, 132,123), and (e, 132,132). 

2 . Always select the left edge. 

Consider (e, e, e). Then {xi, X 2 , Xs)Fo = {x 2 , X 3 , xi+x^) and {xi,X 2 , X 3 )Fi = (xi, X 2 , Xi+ 
X3). Since at the zeroth step we start with {xi,X2,X3) = (1,1,1), it is clear that the 
leftmost triplet will contain the maximal term at each level n since this will lead to 
the greatest rates of growth for each of xi„, X 2 „ and x^^ as n increases. 

Similar arguments can be made for (e, e, 12), (e, e, 23), (e, e, 123), (e, e, 132), (e, 12, e), 

(e, 12,12), (e, 12,13), (e, 12, 23), (e, 12,123), (e, 12,132), and (e, 132, 23). 

3. Alternate between hrst selecting the left edge and then the right edge. 

Consider (e, 13,12). Then (xi, X 2 , X 3 )Fo = (xi+xs, X 3 , X 2 ) and (xi, X 2 , X 3 )Fi = (x 2 , xi, xi + 
X 3 ). Since at the zeroth step we start with (xi,X 2 ,X 3 ) = (1,1,1), it is clear that the 
triplet containing the maximal term at each level will lie on the nodes of the path 
generated by alternating between hrst selecting the left edge and then the right edge. 

Similar arguments can be made for (e, 23, 23), and (e, 123, e). 

□ 

In the above theorem, we have shown that the maximal TRIP-Stern sequence lies on one 
of three possible paths for 26 TRIP-Stern sequences generated by 26 triangle partition maps. 
Using Lemma 9 brings this total up to 26 ■ 6 = 156 maps. 

Question 11. What are the paths for finding maximal terms for the remaining 60 TRIP- 
Stern seguences? 


6.1 Explicit formulas and recurrence relations for sequences of 
maximal terms 


In the above subsection we addressed the positions of maximal terms. Here we present 
formulas and recurrence relations that may be used to hnd the actual values of the maximal 
terms for 120 TRIP-Stern sequences. 


Theorem 12. The maximal term rUn in the TRIP-Stern seguences corresponding to the 
permutation triplets (e, e, 13), (e, 12, e), (e, 12,12), (e, 12,13), (e, 12, 23), (e, 12,123), (e, 12,132), 
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(e, 13,13), (e, 23,13), (e, 123,13), and (e, 132,13) is given by the formula 


m„. = 


-■((75-3) (1 - 75)” + (3 + 75) (1 + 75)”) 

T! 


which corresponds to the Fibonacci recurrence relation = m„_i + m „_2 (AOOOOfS). 

Proof. By Theorem 10, we know that the third term in each of the triplets given by following 
the left-most path in the tree generated by select permutation triplets is a maximal term; 
similarly, the hrst term in each of the triplets given by following the right-most path in 
the tree generated by select permutation triplets is a maximal term. Hence, for sequences 
of maxima found by following the left-most path, all that remains to hnd is to hnd the 
third term in the triplet (1,1, ^)Fo', for each of the permutation triplets listed in the theorem 
that are generated by choosing the left-most path, this third term corresponds to the desired 
explicit formula. Similarly, for sequences of maxima found by following the right-most path, 
all that remains to hnd is to hnd the hrst term in the triplet (1,1, for each of 

the permutation triplets listed in the theorem that are generated by choosing the right-most 
path, this hrst term corresponds to the desired explicit formula. It is easy to check using 
standard methods, as in Matthews [40], that = m„_i -|- m„_ 2 . 

As an example, let us consider (e, e, 13), for which 


Fi = 



Then 


Ff = 


! 2—1 (-(1-^/5)"+' + (l+^/5)"+') 

75 

0 

2-(-(l-^/5)"+(l+^/5)") 

\ v/5 


0 

1 


2-(-(l-V5)" + (l+^)") 

75 

0 

2—1 ((1-^/5)" (l+y5) + (-l+V5) (1+75)") 
75 

SO that 

. /2—l((^/5-3)(l-V5)" + (3+^/5)(l+^/5)") 2—i ((l+^)"+'-(l-75)"+') 

U, -L, ^5 , -L, ^5 

Hence, we can see that 

2—1 (^(VS - 3) (1 - Vs)” + (3 + x/5) (1 + Vs)”) 


\ 


/ 


= 


^/5 


Lastly, it is clear that -|- m „_2 since 


m„_i -h m„_2 = 
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□ 


2-i-(-i)((l-x/5)" '(-3+^/5) + (l+^/5)" ^(3+x/5)) ^ '(-3+^/5) + (l+^/5)" "(3+^/5)) 

Ts ^ 75 

2-"-' ((^5 - 3) (1 - V5)" + (3 + x/5) (1 + \/5)”) 




= 


Theorem 13. The maximal term mn in the TRIP-Stern sequences corresponding to the 
permutation triplets (e, e, e), (e, e, 12), (e, e, 23), (e, e, 123), (e, e, 132), (e, 13,123), (e, 23,123), 
(e, 123,123), and (e, 132,123) is given by the formula 

mn = aiP'f + 02/^3 + “ 3/52 ) 

where the ai’s are roots of 31x^ — 31x^ — 12x — l = 0 while the fdi’s are the roots of x^—x"^ — 1 = 
0. This corresponds to the well-known recurrence relation mn = mn-i + mns fA000930) . 

Proof. The proof, except for some technical details, is identical to the previous proof. By 
Theorem 10, we know that the third term in each of the triplets given by following the left¬ 
most path in the tree generated by select permutation triplets is a maximal term; similarly, 
the first term in each of the triplets given by following the right-most path in the tree 
generated by select permutation triplets is a maximal term. Hence, for sequences of maxima 
found by following the left-most path, all that remains to hnd mn is to find the third term 
in the triplet (1,1, 1)Fq] for each of the permutation triplets listed in the theorem that 
are generated by choosing the left-most path, this third term corresponds to the desired 
explicit formula. Similarly, for sequences of maxima found by following the right-most path, 
all that remains to find is to find the first term in the triplet (1,1,1)7^7; for each of 
the permutation triplets listed in the theorem that are generated by choosing the right-most 
path, this hrst term corresponds to the desired explicit formula. Using standard methods (as 
in Matthews [40]) identical to those used in the example found in the proof of the previous 
theorem, it is easy to see that = mn-i + mns- 

□ 

In the above two theorems, we have presented explicit formulas and recurrence relations 
for the maximal terms of TRIP-Stern sequences generated by 20 maps. Using Lemma 9 
brings this total up to 20 • 6 = 120 maps. 

Question 14. What are the recurrence relations for the other 96 TRIP-Stern sequences? 

Conjecture 15. The following table presents some numerical results and conjectured re¬ 
currence relations for TRIP-Stern maximal terms corresponding to the above-mentioned 96 
maps. 
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First 11 maximal terms 

Conjectured recurrence relation = 

(e,To,ri) 

A-number 

1,2,3,5,7,11,16,25,36,56,81 

unknown 

(e,13,e),(e,123,12) 

A271485 

1,2,3,4,6,9,13,19,28,41,60 

rrin-i -i- m„3 

(e,13,12) 

A000930 

1,2,3,4,6,8,11,16,22,30,43 

unknown 

(e, 13,23), (e, 23,12) 

A271486 

1,2,3,4,6,8,11,17, 23,32,48 

unknown 

(e, 13,132), (e, 132,12) 

A271487 

1,2,3,4,6,8,11,15,21,30,41 

unknown 

(e,23,e),(e,123,23) 

A271488 

1,2,3,4,5,7,9,12,16,21 

mn-2 + rrin-s 

(e, 23,23), (e, 23,132) 

(e, 132, 23), (e, 132,132) 

A000931 

1,2,3,5,8,13,21,34,55,89,144 

rrin-i + mn-2 

(e,123,e) 

A000045 

1,2,3,4, 5, 7,10,13,18,25,34 

unknown 

(e, 123,132), (e, 132, e) 

A271489 


Table 2 : Conjectured recurrence relations for maximal terms 


Note that in the above table we only included maps of the form (e, ro,ri); as before, 
Lemma 9 brings the total up to 16 ■ 6 = 96 maps. 


7 Minimal terms and positions thereof 

We now investigate the positions of the minimal terms in each level for various TRIP-Stern 
sequences. This is a bit easier than the analogous investigation for the maximal terms. 

Theorem 16. The minimal terms bn in the TRIP-Stern sequences corresponding to the 
permutation triplets (e, 12,12), (e, 12,123), (e, 12,132), (e, 13,123), (e, 13,132), (e, 23,12), 
(e, 23,123), (e, 23,132), (e, 123,12), (e, 123,123), and (e, 123,132) lie on the left-most path 
in the corresponding TRIP-Stern tree. 

Proof. For the permutation triplets (e, 13,123) and (e, 13,132), the transformation 


{xi,X2,X3) H- (xi,X2,X3)Fo 

flips the positions of X 2 and X 3 in the triplet. Recall that application of Fq corresponds to 
following the left-most child at each node in the TRIP-Stern tree corresponding to some 
permutation triplet. Therefore, as we start with {xi,X 2 ,X 3 ) = (1,1,1), it is clear that the 
minimal terms of the TRIP-Stern sequences corresponding to these permutation triplets he 
on the left-most path in the corresponding TRIP-Stern tree. 

For the rest of the above permutation triplets, the transformation 


{Xi,X2,X3) H- {Xi,X2,X3)Fo 


leaves either xi, X 2 , or X 3 hxed. As in the hrst part of this proof, since we start with 
{xi,X 2 ,X 3 ) = (1,1,1), it is clear that the minimal terms of the TRIP-Stern sequences cor¬ 
responding to the above permutation triplets lie on the left-most path in the corresponding 
TRIP-Stern tree. 

□ 
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Theorem 17. The minimal terms bn in the TRIP-Stern sequences corresponding to the 
permutation triplets (e, e,e), (e, e, 12), (e,e, 13), (e,e, 23), (e, 13,e), (e, 13,13), (e, 13,23), 
(e, 132,e), (e, 132,12), (e, 132,13), and (e, 132,23) lie on the right-most path in the corre¬ 
sponding TRIP-Stern tree. 

Proof. For the permutation triplets (e, e, 12) and (e, 132,12), the transformation 

{xi,X2,X3) H- {xi,X2,X3)Fi 

flips the positions of xi and X 2 in the triplet. Recall that application of Fi corresponds to 
following the right-most child at each node in the TRIP-Stern tree corresponding to some 
permutation triplet. Therefore, as we start with {xi,X 2 ,X 3 ) = (1,1,1), it is clear that the 
minimal terms of the TRIP-Stern sequences corresponding to the above permutation triplets 
he on the right-most path in the corresponding TRIP-Stern tree. 

For the rest of the above permutation triplets, the transformation 

{Xi,X2,X3) H- {Xi,X2,X3)Fi 

leaves either xi, X 2 , or ^3 hxed. As in the hrst part of this proof, since we start with 
{xi,X 2 ,X 3 ) = (1,1,1), it is clear that the minimal terms of the TRIP-Stern sequences corre¬ 
sponding to the above permutation triplets he on the right-most path in the corresponding 
TRIP-Stern tree. 

□ 

Theorem 18. The minimal terms bn in the TRIP-Stern sequences corresponding to the per¬ 
mutation triplets (e, 12,e), (e, 12,13), (e, 12,23), (e, 13,12), (e, 23,e), (e, 23,13), (e,23,23), 
(e, 123, e), (e, 123,13), and (e, 123,23) lie on both the right-most and left-most paths in the 
corresponding TRIP-Stern tree. 

Proof. For the permutation triplet (e, 13,12), the transformation (xi, X 2 , X 3 ) (xi, X 2 , X 3 )Fq 
flips the positions of X 2 and X 3 in the triplet and the transformation (xi, X 2 , X 3 ) i-4 (xi, X 2 , X 3 )Fi 
flips the positions of xi and X 2 . As we start with (xi,X 2 ,X 3 ) = ( 1 , 1 , 1 ), it is clear that the 
minimal terms of the TRIP-Stern sequences corresponding to the above permutation triplets 
lie on both the left- and right-most paths in the corresponding TRIP-Stern tree. 

For the rest of the above permutation triplets, both the transformations (xi,X 2 ,X 3 ) 1 -^ 
(xi,X 2 ,X 3 )Fo and (xi,X 2 ,X 3 ) 1 -^ (xi,X 2 ,X 3 )Fi leave either xi, X 2 , or X 3 hxed. Since we start 
with (xi,X 2 ,X 3 ) = (1,1,1), it is clear that the minimal terms of the TRIP-Stern sequences 
corresponding to the above permutation triplets lie on the right-most and left-most paths in 
the corresponding TRIP-Stern tree. □ 

Corollary 19. The minimal terms bn in the TRIP-Stern sequences corresponding to the 
permutation triplets mentioned in Theorems 16, 17, and 18 all have value 1. 

Proof. This follows immediately since in each case we start with (xi,X 2 ,X 3 ) = (1,1,1), and 
at least one of the components of this triplet gets carried over to at least one triplet in the 
next level of the corresponding TRIP-Stern sequence by the action of Fq or Fi, as outlined 
in the proofs of the theorems mentioned in the corollary. □ 
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In the above theorems, we have found the values and positions of the minimal terms of 
TRIP-Stern sequences generated by 32 maps. The results of Lemma 9 bring this total up to 
32 • 6 = 192 maps. 

Question 20. What are the values and positions of the minimal terms of TRIP-Stern se¬ 
quences generated by the remaining 2f maps? 

8 Level sums 

The following section examines the sums of the entries in each level of a TRIP-Stern sequence, 
in direct analogue to the level sums found by Stern for his diatomic sequence (Lehmer [38] 
presents this as Property 2). As the level n grows large, the ratio between the sums of the 
entries in successive entries approaches an algebraic number of degree at most 3. The ratios 
between the hrst, second, and third coordinates of a given level approach ratios in the same 
number held. This section provides a closed form for the sums of the entries in each level. 

Definition 21. Consider the TRIP-Stern sequence for arbitrary [a, tq, ti). Let Si{n) be the 
sum of the hrst entries of the triples in the level, let 5*2 (n) be the sum of the second 
entries, and let S' 3 (n) be the sum of the third entries. The sum of all entries in a given level 
is S{n) = S'i(n) 5*2 (n) S' 3 (n). 

Proposition 22. For each n, 

{Si{n), S 2 {n), Ss{n)) = {S,{n - 1), S^in - 1), Ss{n - 1)) (Fq + F,). 

Proof. The base case follows by dehnition. Say that we are at the level. In order to 
generate the next level of triplets, we apply Fq and Fi to each triplet in the level. It is 
clear that Si{n + 1) is obtained by taking the sum of the hrst components of OmF^ + a^Fi 
over all triplets in the level; similarly for S' 2 (n -|- 1) and S^i^n -\- 1). As a result, it is 
clear that 

(^i(n), S 2 {n), Ssin)) = {S,{n - 1), S 2 {n - 1), S^in - 1)) (Fq + F,). 

□ 


8.1 Level sums for (e, e, e) 

Proposition 23. Let a be the real zero of + 5x — 4. Then 


for 1 < i < 3 and 


1- Si{n) 

™ Si{n - 1) 


= a 


lim 

n^oo 


S{n) 
S{n - 1) 


= a. 
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Proof. The characteristic polynomial of Aq + Ai is + 5x — 4. This polynomial has 

real zero a ~ 2.69562, and two complex zeros of smaller absolute value. Thus, as n —)> oo, 
the vector 

(^i(n), = (1,1, l)(dlo + Ai)" 

approaches the eigenvector a corresponding to the eigenvalue a. So as n approaches infinity, 
the vector (1,1, l)(Ao + Ai)^ approaches the subspace generated by a. Hence, 

r Si{n) 

n^oo bi{n — 1) 


for 1 < i < 3 and 


as claimed. 


lim 

n^oo 


S{n) 
S{n - 1) 


= a, 


□ 


Proposition 24. VTe have 

,. S2{n) S^in) 

hm =0 — 1 and hm = a — 1. 

n^oo Si{n) n^oo 62 

Proof. The actions of Aq and Ai yield the recurrence relation S'i(n + 1) = S'i(n) + 5*2 (n). 
Dividing by S'i(n) gives 

>5'i(n + 1) _ , »S'2(n) 

S'i(n) Si{n) 

Taking the limit as n —)■ oo and applying Proposition 23 yields 

= 0 — 1 . 

1 ) = S 2 {n) + Sq^h) yields 
= 0 — 1 . 

□ 


S2{n) 


lim . . 

n^oo ibi[n) 


Similarly, using the recurrence relation 5*2 (n + 

Ssin) 


lim . . 

n^oo D2[n) 


8.2 Level sums for arbitrary (cr, tq, ri) 

The properties of level sums for an arbitrary TRIP-Stern sequence are similar to those of the 
TRIP-Stern sequence for (e,e, e). For any (a, ro,ri), it can be shown by computation that 
Fq + Fi has an eigenvalue of largest absolute value. In fact, this eigenvalue is an element of 
the interval [2, 3]. 
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Proposition 25. For any TRIP-Stern sequence, we have 


lim 

n—^oo 


Sjjn) 

Si{n - 1 ) 


= a 


for 1 < i < 3 and 


lim 


S{n) 


= a, 


S{n - 1) 

where a is the eigenvalue of Fq + Fi of largest absolute value. Furthermore, a is an algebraic 
number of degree at most three, and a G [2, 3]. 


Analogous recurrence relations give relations between the limits of Si{n), S 2 {n), and 
Ssi^n). These are not always as clean as in the case of the TRIP-Stern sequence for (e, e, e), 
but the ratios between these limits are contained in the number held Q{a). 

We next examine the closed forms for S{n). 


Theorem 26. The family of triangle partition maps leads to 11 distinct sequences of sums 
with recurrence relations and explicit forms as shown in the tables below. 

Proof. The proof follows by direct calculation. For each triangle partition map T„^ro,Tn 
compute a modihed TRIP-stern sequence given by setting oi = (a, b, c) instead of setting 
oi = (1,1,1) as we had done in Section 4, partitioning the sequence into levels as before. 
Sum the terms of each level n to yield a sequence of row sums (5'(n))^>p 

If we can prove that the hrst m terms of satisfy an (m — l)-term recurrence 

relation, it follows that the sequence must be generated by that recurrence relation. We 
have carried out this procedure for all 216 permutation triplets (a, ro,ri) to hnd recurrence 
relations for the associated row sums, from which the explicit form for the term in the 
sequence (*S'(?7,))^>^ was easily calculated. The results are presented in the tables below - 
indeed, the family of triangle partition maps generates only 11 distinct row sums. The hrst 
column lists the recurrence relation, the second lists the explicit form of that recurrence 
relation, and the third lists the permutation triplets whose TRIP-Stern level sums follow 
this relation. Note that Greek letters represent zeros of certain polynomials; see the key 
below. For example. Root [29x^ — 87x^ — 5x — 1,1] —)■ means “let ai be the hrst root of 
29x3 - 87x2 _ 53. _ 1 = 0.” □ 




Recurrence relation for S{n) 

Explicit form for Sin) 

(e,To,ri) 

4S{n — 3) — 55'(n — 2) + 4S{n — 1) 

+ CX 2 P 2 4~ CXsPs 

(e,e,e), 

(e, 123,123) 

25(n-2) + 25(n-l) 

|((9-5V3) (l-\/3)" 

+ (1 + ^3) (9 + 5^)) 

(e,e, 12), 

(e,e, 123), 
(e,13,12), 

(e, 13,123) 

5(n-3)-5(n-2) + 35(n-l) 

+ 73^2 + 72^3 

(e,e, 13), 

(e, 12,123) 

-S{n - 3) + 2S{n - 2) + 2S{n - 1) 

2-'*(-2(3-75)”(-2+75) + (3+75)"(i1+575)) 

5+75 

(e,e,23), 

(e,12,23), 

(e, 12,132), 
(e,23,e), 
(e,23,13), 

(e, 23,123), 

(e, 123,23), 

(e, 123,132), 

(e, 132, e), 

(e, 132,13) 

6S{n — 3) + 2S{n — 2) + S{n — 1) 

eiCr + e2C? + e3C? 

(e,e, 132), 

(e, 132,123) 

5S{n — 1) — 65'(n — 2) 

2” + 2 3” 

(e, 12, e), 
(e,12,13), 

(e, 123,e), 

(e, 123,13), 

—4S{n — 3) + S{n — 2) + 3S{n — 1) 

rjiO^ + r?20^ + 73^3 

(e,12,12), 

(e,13,13) 

-S{n - 3) - 3S{n - 2) + 4S{n - 1) 

ilKi + ^2^2 T ^3^3 

(e, 13, e) 

(e, 123,12) 

-6S{n - 3) + 4S{n - 2) + 2S{n - 1) 

-^ 1/^1 + ^ 2/^2 + 

(e,13,23), 

(e,23,12) 

5(n-3) + 45(n-2) + 5(n-l) 

yiQ + ^2^2 + 

(e, 13,132), 

(e, 132,12) 

4S{n - 2) + Sin - 1) 

i((51-1371^(^1-^)" 

+ (i(l + ^) (51 + 13^) 

(e,23,23), 

(e, 23,132), 

(e, 132, 23), 

(e,132,132) 


Table 3: Level sums 
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Root [29x^ — 87x'^ — 5x — 1,1] —^ cti 
Root [29x^ — 87x^ — 5x — 1, 3] —^ as 
Root [x^ — 4x^ + 5x — 4, 2] —>• /32 
Root [76x^ — 228x^ + 28x — 1,1] —>■ 71 
Root [76x^ — 228x^ + 28x — 1, 3] —>• 73 
Root [x^ — 3x^ + X — 1, 2] —)■ ^2 
Root [1176x^ — 3528x^ — 119x — 1,1] —^ ei 
Root [1176x^ — 3528x^ — 119x —1,3]—^ 63 
Root [x^ — x^ — 2x — 6,2] —(2 
Root [229x^ — 687x^ + 230x + 4,1] —rji 
Root [229x^ — 687x^ + 230x + 4, 3] —>■ 773 
Root [x^ — 3x^ — X + 4, 2] —> 02 
Root [49x^ — 147x^ + 35x — 1,1] —ii 
Root [49x^ — 147x^ + 35x — 1, 3] —>■ 7.3 
Root [x^ — 4x^ + 3x + 1, 2] —>■ K2 
Root [404x3 - 1212x2 + 24x + 1,1]^ Ai 
Root [404x3 — 1212x2 + 24x + 1, 3] —>■ A3 
Root [x3 — 2x2 — 4x + 6, 2] —^ 7x2 
Root [169x3 — 507x2 ^ 

Root [169x3 — 507x2 ^ 3 j 

Root [x3 — x2 — 4x — 1, 2] —^ ,^2 


Root [29x3 — 87x2 _ 5^; _ 2] —>■ a2 

Root [x3 — 4x2 + 5x — 4,1] —>• f3i 
Root [x3 — 4x2 + 5x — 4,3] — f3s 
Root [76x3 _ 228x2 + 28x — 1,2]—^ 72 
Root [x3 — 3x2 _|_ 2- _ 1] 

Root [x3 — 3x2 _|_ 2- _ X^ 3] 

Root [1176x3 — 3528x2 — 119x — 1,2] —> €2 

Root [x3 — x2 — 2x — 6,1] —>• Cl 

Root [x3 — x2 — 2x — 6,3] —>■ Cs 

Root [229x3 — 687x2 _|_ 230x + 4, 2] —>■ r/2 

Root [x3 — 3x2 _ 2; _j_ 4^ X] —>• 01 

Root [x3 — 3x2 _ 2- _j_ 4^ 3j 

Root [49x3 — 147x2 + 35x — 1, 2] —>• 12 
Root [x3 — 4x2 + 3x + 1,1] —)• Ki 
Root [x3 — 4x2 + 3x + 1, 3] —)• K 3 
Root [404x3 — 1212x2 + 24x 4- 1,2]—^ A 2 
Root [x3 — 2x2 — 4x + 6 ,1] —>■ /ii 
Root [x3 — 2x2 — 4x + 6 , 3] —>■ 

Root [169x3 — 507x2 _|_ x^ 2] —>• 1^2 
Root [x3 — x2 — 4x — 1,1] —!> Cl 
Root [x3 — x2 — 4x — 1, 3] —^ Cs 


Table 4: Key to level sums table 


Note that in the above table we only included maps of the form (e, ro,ri); as before, 
Lemma 9 brings the total up to 36 * 6 = 216 maps. 


8.2.1 Six row sum recurrence relations already well-known 

It turns out that 6 of the 11 unique row sum recurrence relations or sequences - either 
with the same initial terms or with different ones - had already been placed on Sloane’s 
encyclopedia; in particular, these recurrence relations and corresponding Sloane sequence 
numbers are as follows: 

S{n) = 2S{n - 1) + 2S{n - 2) f 21080040 . 21155020 1 

^(n) = 3^(n - 1) - S{n - 2) + S{n - 3) 1 21200752 1 

S{n) = 2S{n-l)+ 2S{n-2) - S{n-3) 1 21061646 1 

S{n) = 5S{n - 1) - 6S{n - 2) 1 21007689 1 

^(n) = 4^(n - 1) - 3^(n - 2) - - 3) 1 21215404 1 

S{n) = S{n — 1) + 4:S{n — 2) 1 A006131 1 
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The new sequences^ are generated by the following recurrence relations: 


S{n) = 4S'(n — 1) — 5S'(n — 2) + 4S'(n — 3) (74278612) 

S{n) = S{n-1)+ 2S{n-2)+ 6S{n-3) 1 74278613 1 
S{n) = 3S{n-l) + S{n-2)-4S{n-3) (A2786U ) 
S{n) = 2S{n - 1) + 4^(n - 2) - 6^(n - 3) 1 74278615 1 

S{n) = S'(n — 1) + 4S'(n — 2) + S'(n — 3) 1 74278616 1 


9 Forbidden triples for (e, e, e) 

This section explores which points in 7? appear as terms in the TRIP-Stern sequence for 
(e, e, e) and which points do not. We call these latter points forbidden triples. 

Definition 27. Let S denote the set of points given by the TRIP-Stern sequence for (e, e, e), 
let P = {{x,y, z) & < X < y < z} U (1,1,1) denote the set of potential entries in S and 

let F = P \ S' denote the set of forbidden triples. 

Proposition 28. We have that S P. 

Proof. Suppose {a,b,c) G S. We proceed by induction on the level of {a,b,c). If {a,b,c) 
is in the second level, then {a,b,c) = (1,1,2), so it satisfies 0 < a < 6 < c. Suppose 
that 0 < a < fe < c for all elements {a,b,c) of the level. If {a,b,c) is in level n -|- 1, 
then (a, b, c) = (d, e, f)Ai^ with in G {0,1} for some (d, e, f) satisfying 0 < d < e < f. If 
in = 0, then (a, b, c) = (e, f,d + f). We know that 0 < e < /, and f < d + f because d is 
positive. If in = 1, then (a, b, c) = (d, e,d+ f). By the inductive hypothesis, 0 < d < e and 
e<f<d + f. □ 

By Proposition 28, elements of F are precisely the potential entries that are not in the 
TRIP-Stern sequence for (e,e, e). 

Definition 29. We define the inverse map G of a triple {a,b,c) to be 


G 


{a,b,c — a), ii a + b < c; 

{c — b,a,b), ii a + b > c and a < b, or a = 1 = c — b. 


The map G is 74^^ if a -|- 5 < c and is 74 q ^ if a -|- 6 > c is 74 q However, G is not defined 
on all points of P. We will show that, for any point (a, b, c) in P that is not (1,1, 2), at most 
one of (a, b, c)Af^ or (a, b, c)Aq^ can lie in the TRIP-Stern sequence for (e, e, e). 

Proposition 30. If {a, b, c) G P, and (a, b, c) ^ (1,1, 2), then either (a, b, c)Aq^ or (a, b, c)Af^ 
is not in P. 

^We have placed these sequences on Sloane’s encyclopedia. 
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Proof. First assume (a, b, c) has a + b > c. Then 

(a, b, c)Af^ = (a, b,c — a) 

However, we have that b > c — a. Then (a, b, c)Af^ is not in P, unless (a, b,c — a) = (1,1,1), 
which occurs if and only if a = 6 = 1, c = 2. Now suppose (a, b, c) has a + b < c. Then 

(a, b, c)Aq^ = {c — b, a, b). 

In order for this to lie in P we need c — b < a < b. In particular, c < b + a, contradicting 
the assumption. Thus, either (a, b, c) = (1,1, 2) or at most one of (a, b, c)Af^ or (a, b, c)Hq ^ 
lies in P. □ 

Corollary 31. For every X E S except (1,1,1), the point X appears exactly twice in the 
TRIP-Stern sequence for (e,e, e). Furthermore, G maps X to the unique Y E S such that 
either Y Aq = X or YAi = X. 

Proof. If a point (a, b, c) ^ (1,1,1) lies in S', then by dehnition either (a, b, c)Af^ or (a, b, c)Aq^ 
must lie in S. Proposition 30 implies that only one of these can be in P. Thus, exactly one 
of these points is in S. This takes care of the second statement. 

Because the left and right subtrees of the TRIP-Stern sequence for (e, e, e) are symmetric, 
we need only show that each X E S appears exactly once in the set of points generated by 
the action of Ai and Aq on (1,1,2). Now suppose that up to level n each entry in the 
TRIP-Stern sequence appears only once. Then, in level n -|- 1, each element X corresponds 
to either XAf^ or XAq^ in level n. By Proposition 30, exactly one of XAf^ or XAq^ will 
he in level n, and this is the only element that goes to X under one of Aq or Ai. This is 
the unique preimage in S, under Aq and Hi, that goes to X. Then each entry in level n -|- 1 
appears for the first time, and in the level n -|- 1 there are no repeated entries, completing 
the induction. □ 

Definition 32. We define a germ to be any element {a,a,b) E P with b < 2a. We call 
the set of all elements generated by action of Aq and Hi on (a, a, b) the tree generated by 
{a,a,b). 

Observe that the tree generated by (1,1,1) is precisely the TRIP-Stern Sequence for 
(e, e, e). The only points for which G is not defined in P are germs. Moreover, each applica¬ 
tion of G to X E P decreases (strictly) the sum of the entries of X. Since G is well-defined on 
all of P, excluding germs, after some number of applications of G to an element X, we find 
a germ that generates X. The following lemma, whose proof is straightforward, is needed to 
strengthen the result. 

Lemma 33. For all X E P, we have G((X)Hi) = X and G((X)Ho) = X. 

This lemma shows that there is a unique germ generating X for each X E P. 
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Definition 34. Let the germ of X be the value G^{X) such that n is the largest integer 
for which G'^{X) is defined. As noted above, this G^{X) will be a germ. 

Theorem 35. Every element of P lies in exactly one tree generated by a germ. Furthermore, 
an element X E P lies in the tree generated by Xq if and only if the germ of X is Xq. In 
particular, one can determine the germ of any given triple (a, b, c) in a finite number of steps. 

Proof. If the germ of X is Xq, then since G acts as or Aq^ at each step, we have that 
X can be written as (Xo)Aj^Aj 2 • • • Ai^ for some ij G {0,1}. Conversely, if X is in the tree 
generated by Xq, then we can write X as (A'o)AijAj2 • • • for some ij G {0,1}. But then 
by Lemma 33, we have that G”(X) = Xq, so the germ of X is Xq. □ 

Corollary 36. No germ except (1,1,1) lies in the TRIP-Stern seguence for (e,e,e). All 
elements that do not lie in the TRIP-Stern Seguence are given by the action of Aq and Ai 
on a germ other than (1,1,1). 

Proof. Both Af^ and Aq ^ take elements of the form (a, a, b) with b < 2a outside of P. Then 
(a, a, b) cannot be reached by the action of Ai or Aq on an element of P. □ 

No doubt a similar analysis can be done for any TRIP-Stern sequence, though we do not 
know how clean the analogues would be. 


10 Generalized TRIP-Stern sequences 

In the original definition of a TRIP-Stern sequence from Section 4.1, we set the first triple 
in the sequence to be oi = (1,1,1). However, there is nothing canonical about this choice, 
which leads us to construct generalized TRIP-Stern sequences, where we set the initial triple 
to be oi = (a, b, c), for some a,b,cE M. 

10.1 Construction of generalized TRIP-Stern sequences 

Definition 37. For any permutation triplet (a, ro,ri) in Af, the generalized TRIP-Stern 
seguence of (a, ro,ri) is the unique sequence such that, for some a,b,c G M, Oi = {a,b,c) 
and, for n > 1, 

r ®2n ■ Pq 

{ fl2n+l = 0.n ' Fl 

The level of the generalized TRIP-Stern sequence is the set of with 2"'“^ < m < 2”. 

As for the standard TRIP-Stern sequence, each choice of (a, ro,ri) also produces some 
generalized TRIP-Stern sequence. 

Definition 38. Let {a,b,c) be any triple of real numbers and let (a, ro,ri) E Sq x S 3 x S 3 . 
Let T{a,TQ,Ti)(^a,b,c) denote the tree generated from (a, ro,ri) using {a,b,c) as a seed. 
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10.2 Maximum terms and positions thereof for generalized TRIP- 
Stern sequences 

This section examines maximum terms in any given level of a generalized TRIP-Stern se¬ 
quence, as well as the positions of those maximum terms within the given level. For a given 
seed (a, b, c) and for an n-tuple v of zeros and ones, dehne 

A{v) = 

which can be written as A(n) = (&i,& 2 ,^ 3 )- As before, let |n| denote the number of entries 

in V. 

Definition 39. The maximum entry on level n of a generalized TRIP-Stern sequence is 
= max|^l=„ maxi6{i,2,3} 

Lemma 40. Suppose (a, Tq, Ti) G Rf, k G S 3 , and let v = (ii,..., A). Then 

1) (^) c) ' • CxAjjTjj ' ' ' ^ 

Proof. The proof uses the same technique as in Lemma 9. □ 

We will now examine the sequences of generalized TRIP-Stern maximal terms induced 
by select triangle partition maps. 

We hnd two broad classes of generalized TRIP-Stern sequences for paths through the 
corresponding trees that will locate the maximal terms at each level. These depend to some 
extent on the initial seeds. As the proofs are straightforward and similar to the earlier ones, 
we will omit them. 

Theorem 41. We have 

1. For any seed {a,b,c), if a > b > c > 0, the sequence of maximal terms for TRIP- 
Stern sequences induced by the nine maps (e, 13,123), (e, e, 13), (e, 13,13), (e, 23,13), 
(e, 23,123), (e, 123,13), (e, 123,123), (e, 132,13), and (e, 132,123) lies on the con¬ 
nected path of the tree T{e,To,Ti)(a,b,c) rnade by always selecting the right edge of 

T (e, Tq, Ti) (^a,h,c) ■ 

2. For any seed {a,b,c), if 0 < a < b < c, the sequence of maximal terms for TRIP-Stern 
sequences induced by the eleven maps (e, e, e), (e, e, 12), (e, e, 23), (e, e, 123), (e, e, 132), 
(e, 12,e), (e, 12,12), (e, 12,13), (e, 12,23), (e, 12,123) and (e, 12,132) lies on the con¬ 
nected path through the tree T{e,To,Ti)(^a,b,c) made by always selecting the left edge of 

Ri R)(a,b,c) ■ 

In the above theorem, under select conditions, we have accounted for sequences of gen¬ 
eralized TRIP-Stern sequence maximal terms generated by 20 maps. Lemma 40 brings this 
total up to 20 ■ 6 = 120 maps, as long as the conditions - which simply guarantee that the 
components of the initial seed will have the right magnitudes to satisfy the conditions of 
Theorem 41 after being acted upon by the first k in Lemma 40 - listed in the theorem below 
are satished: 
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Theorem 42. 1. Let {a, Tq, Ti) G S'! be one of the 9 permutation triplets listed in Theorem 

fl.l. For any a seed (a, b, c), the sequence of maximal terms for TRIP-Stern sequences 
induced by maps of the form (kct, roK“^, lies on the connected path through the 

tree roK“^, riK“^)(a,6,c) niade by selecting the right edge ofT^na, TiK~^)(a,b,c), 

as long as the following conditions are satisfied: 


(a) For K = 

(b) for K = 

(c) for K = 

(d) for K = 

(e) for n = 


12, require b > a > c> 0, 

13, require c > fe > a > 0, 

23, require a > c > b > 0, 

123, require c > a > b > 0, and 
132, require h > c > a > 0. 


2. Let (a, ro,ri) G S'! be one of the 11 permutation triplets listed in Theorem fl.2. For 
any seed {a,b,c), the sequence of maximal terms for TRIP-Stern sequences induced 
by maps of the form {Ka,TQK~^,TiK~^) lies on the connected path through the tree 
,TiK~^)f^a,b,c) made by selecting the left edge of Tinner, ToK,~^,TiK~^)(^a,b,c), 
long as the following conditions are satisfied: 


(a) For K = 

(b) for n = 

(c) for K = 

(d) for K = 

(e) for K = 


12, require 0 < b < a < c, 

13, require 0 < c < 6 < a, 

23, require 0 < a < c < 6, 

123, require 0 < c < a < b, and 
132, require 0 < b < c < a. 


Proof. This follows immediately by the results of Theorem 41 and remembering that k, is 
simply a permutation. □ 


For any node (r, s,t) in the tree T(a, tq, ri)(„,&,c), it is natural to consider the positions 
of maximal terms within the tree that would be generated using that (r, s,t), and not the 
original (a, b, c) as its seed. Clearly, the problem of characterizing these terms is equivalent 
to characterizing the maximal terms on the nodes below and connected with (r, s, t); in this 
sense, the problem is one of hnding a sequence of local maximal terms. It is straightforward 
to prove the analogous theorems. 


10.3 Minimal terms and positions thereof 

We also have analogs for hnding minimal terms for a number of generalized TRIP-Stern 
sequences. As the proofs are similar to the earlier ones, we omit them. 

Theorem 43. For any seed {a,b,c), the minimal terms bn in the TRIP-Stern sequences 
corresponding to the permutation triplets listed below lie on the left-most path in the corre¬ 
sponding TRIP-Stern tree: 
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1. For the maps (e, 12, e), (e, 12,12), (e, 12,13), (e, 12, 23), (e, 12,123), (e, 12,132) under 
the condition that h < a and b < c - the minimal term will have value b at every 
level. 

2. For the maps (e,123,e), (e, 123,12), (e, 123,13), (e, 123, 23), (e, 123,123), (e, 123,132) 
under the condition that b < a and b < c or c < a and c < b; correspondingly, the 
minimal term will have value b or c at every level. 

3. For the maps (e, 23, e), (e, 23,12), (e, 23,13), (e, 23, 23), (e, 23,123), (e, 23,132), under 
the condition that c < a and c < b - the minimal term will have value c at every level. 

Theorem 44. For any seed {a,b,c), the minimal terms bn in the TRIP-Stern sequences 
corresponding to the permutation triplets listed below lie on the right-most path in the corre¬ 
sponding TRIP-Stern tree: 

1. For the maps (e, e, 23), (e, 12, 23), (e, 13, 23), (e, 23, 23), (e, 123, 23), (e, 132, 23), under 
the condition that a < b and a < c - the minimal term will have value a at every level. 

2. For the maps (e, e, 13), (e, 12,13), (e, 13,13), (e, 23,13), (e, 123,13), (e, 132,13), under 
the condition that b < a and b < c - the minimal term will have value b at every level. 

3. For the maps (e, e, e), (e, 12, e), (e, 13, e), (e, 23, e), (e, 123, e), (e, 132, e), under the con¬ 
dition that a < b and a < c XOR b < a and b < c; correspondingly, the minimal term 
will have value a or b at every level. 

In the above two theorems, we have found the positions and values of the minimal terms 
of generalized TRIP-Stern sequences generated by 27 maps under certain conditions on the 
initial seed. Note that certain maps appear in both the “left” and “right” lists. Using the 
results of Lemma 40 and imposing conditions analogous to those used in Theorem 42 for 
maximal terms - which simply guarantee that the components of the initial seed will have 
the right magnitudes to satisfy the conditions of the above two theorems after being acted 
upon by the hrst k in Lemma 40 - brings this total up to 27 • 6 = 162 maps. 

10.4 Level sums for generalized TRIP-Stern sequences 

It is natural, as was done with standard TRIP-Stern sequences, to consider level sums for 
generalized TRIP Stern sequences. We take an identical approach to that for standard 
TRIP-Stern sequences, though the explicit forms are more complex, as one would expect. 

Theorem 45. The family of triangle partition maps leads to 11 distinct sequences of sums 
{S{n))n>i with recurrence relations and explicit forms as shown in the tables below. 
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Proof. The proof follows by direct calculation. For each triangle partition map To-_t-o,ti, 
compute a generalized TRIP-stern sequence given by setting oi = (a, b, c) instead of setting 
Oi = (1,1,1) as we had done in Section 4, partitioning the sequence into levels as before. 
Sum the terms of each level n to yield a sequence of row sums (-S'(n))^>p 

If we can prove that the first m terms of satisfy an (m — l)-term recurrence 

relation, it follows that the sequence must be generated by that recurrence relation. We 
have carried out this procedure for all 216 permutation triplets (a, ro,ri) to hnd recurrence 
relations for the associated row sums, from which the explicit form for the term in the 
sequence (S(n))^^j^ was easily calculated. 

The results are presented in the tables below - indeed, the family of triangle partition 
maps generates only 11 distinct row sums. The hrst column lists the recurrence relation, the 
second lists the explicit form of that recurrence relation and the third lists the permutation 
triplets whose TRIP-Stern level sums follow this relation. Note that Greek letters represent 
zeros of certain polynomials; see the key below. For example. Root — 4x^ + 5x — 4,1] —> 
ai means “let ai be the first root of — Ax"^ + 5x — 4 = 0.” □ 
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Recurrence relation for S{n) 

Explicit form for S{n) 

(e,To,Ti) 

45(n - 3) - 55(n - 2) + 45(n - 1) 

(cCi + brji) a^+a {Pia^ + 

+03 (&73 + 082 ) £2 + 02 (cCa + br] 2 ) 

(e,e,e), 

(e, 123,123) 

2S(n-2)+2S'(n-l) 

i((l-V3)'‘((3-2V3)a 
- (-3 + a/3) 6 + (3 - 2a/3) c) 

+ (1 +V3)”((3 + 2vi)a 

+ (3 + a/s) 6 + (3 + 2a/3) c)) 

(e,e, 12), 
(e,e, 123), 
(e,13,12), 

(e, 13,123) 

S{n — 3) — S{n — 2) + 3S{n — 1) 

aiiOi + bKiOi + cfiiOi + a^2 ^2 
+6^2 ^3 ^3 (^^3 “1“ ^^2 “1^ CM2) C^2M3 

(e,e, 13), 

(e, 12,123) 

-S{n - 3) + 2S{n - 2) + 2S{n - 1) 

2)" (5 +a/ 5) (6 c) 

+ (3-V^)"(-5(-l +V^)a 
+ (5 - 3a/5) & - 2 (-5 + a/5) c) 

+ (3 +V5)"(10 (2 +Vs) a 
+ (15 + 7 a/5) & + 4 (5 + 2V5) c)) 

(e,e,23), 

(e,12,23), 

(e, 12,132), 
(e,23,e), 
(e,23,13), 

(e, 23,123), 

(e, 123,23), 

(e, 123,132), 

(e, 132, e), 

(e, 132,13) 

6S{n - 3) + 2S'(n - 2) + S{n - 1) 

+ boiV^ + apiV^ + 01 ^ 3^2 + bvj}03+ 

1^2 (Vs + b 02 + ap 2 ) + au^pi 

(e,e, 132), 

(e, 132,123) 

5S{n — 1) — QS{n — 2) 

2"& + 3”(a + c) 

(e, 12,e), 
(e,12,13), 

(e, 123, e), 

(e, 123,13) 

-4S’(n - 3) + S{n - 2) + 35(n - 1) 

cr2cr" + avi(7i + bpia^ + ca^ n + acrj £^2+ 
6(72/>2 + 0-3 (ct3 + avs + bps) 

(e,12,12), 

(e,13,13) 

-S{n - 3) - 35(n - 2) + 45(n - 1) 

apiXi + auJ2Xi + 1 X 1 + 0 X 2^2 + o-x^uJi 

+&X2 T2 + Xs (c/’3 + 0^3 + bF 3) 

(e, 13,e), 

(e, 123,12) 

-6S{n - 3) + 45(n - 2) + 25(n - 1) 

a£2^i + ctJiII" + 6<;ini + asilf^ + C'd2^2 

+ ^21^2 + ( 0^3 + C'i?3 + hc,s) 

(e,13,23), 

(e,23,12) 

S{n-3)+iS{n-2) +S{n-1) 

aw2>Ci + cpi>f" + bp2>ti + ax^zui + CK 2 Q 2 
-Fbx^^pi + Xs {azus + cps + bips) 

(e, 13,132), 

(e, 132,12) 

4S{n-2) + S{n-1) 

^((|(l-Vl7))'‘((l7-5Vl7)a 
+ (17 - 3VI7) 6 + (17 - 5VT7) c) 

+ (i (1 + VT7))” ((17 + 5VT7) a 
+ (17 + 3a/T7) 6 + (17 + 5a/T7) c)) 

(e,23,23), 

(e, 23,132), 

(e, 132,23), 

(e, 132,132) 


Table 5: Level sums 
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Root [x^ — 4x^ + 5x — 4,1] —>• «! 

Root [x^ — 4x^ + Sx — 4,3] —)■ as 
Root [582;3 - 58x2 _ -2,2]^ ^2 

Root [x® + 4x2 + X + 2,3] —>■ 73 
Root [116x^ + X + 1,2] ^ £2 
Root [116x3 - 116x2 - 7x - 1,3] ^ Cs 
Root [116x3 - 116x2 + 25x - 2,2] 772 

Root [x3 — 3x2 _|_ 2; — 1^ 2] —>• 02 
Root [19x3 — 19x2 — 3x — 1,1] —>■ ti 
Root [19x3 — 19x2 — 3x — 1,3] —>■ t3 
Root [38x3 _ -(- 10a; — 1, 2] —>• K2 

Root [x3 — 3x2 ^ 2; — 1,1] Ai 

Root [x3 — 3x2 ^ 2; _ 3] A3 

Root [76x3 _ 75^,2 _ 2a; - 1, 2] Ii2 
Root [x3 — x2 — 2x — 6,1] ^ Vi 
Root [x3 — x2 — 2x — 6,3] —>■ 723 
Root [147x3 - 147x2 - 7x - 1,2] ^ ^2 
Root [588x3 _ 5382,2 ^ 773, _ 4^ 1] 

Root [588x3 _ 5882.2 - 4,3] 03 

Root [1176x3 - 1176x2 - 161x - 6,2] ^ p 2 
Root [x3 — 3x2 _ 2; -I- 4^ 1] cr^ 

Root [x3 — 3x2 _ 2; 4- 4^ 3] fjg 

Root [229x3 - 229x2 - 33x + 1,2] T2 
Root [229x3 - 229x2 + 5x + 2,1] ^ 

Root [229x3 — 229x2 + 5x + 2,3] —>■ 773 
Root [229x3 - 229x2 + 61x - 2,2]^ (j >2 
Root [x3 — 4x2 + 3x + 1,1] ^ Xi 
Root [x3 — 4x2 + 3x + 1,3] —>■ Xs 
Root [49x3 — 49x2 + 1,2] —>• 7/^2 
Root [49x3 — 49x2 — 14x + 1,1] —>■ wi 
Root [49x3 — 49x2 — 14x + 1,3] —>■ W3 
Root [49x3 — 49x2 + 14x — 1, 2] —>• Fa 
Root [x3 — 2x2 — 4x + 6,1] ^ Hi 
Root [x3 — 2x2 — 4x + 6,3] Ha 
Root [101x3 — 101x2 + 19X — 1, 2] —>■ <72 

Root [202x3 - 202x2 - 42x - 1,1] ^ ej 

Root [202x3 - 202x2 - 42x - 1, 3] £3 

Root [404x3 — 404x2 — 14x + 3,2] 2 

Root [x3 — x2 — 4x — 1, 1] Ki 
Root [x3 — x2 — 4x — 1,3] —>■ X3 
Root [169x3 — 169x2 — 26x + 1,2] zu 2 
Root [169x3 — 169x2 — 13x + 5,1] —>■ Qi 
Root [169x3 — 169x2 — 13x + 5,3] —?> 

Root [169x3 — 169x2 + 26x — 1,2] —>• 752 


Root [x3 — 4x2 + 5x — 4,2] ^ aa 
Root [58x3 _ 583.2 _ 773 . _ 2 ,1] /3i 

Root [58x3 _ 583-2 _ 773 , _ 2 ,3] ^ ^3 
Root [x3 + 5x2 — 3x + 1,2] —>• S 2 
Root [116x3 - 116x2 - 7x - 1,1] ^ Cl 
Root [116x3 — 116x2 + 25x — 2,1] ^ r]i 
Root [x3 — 3x2 _|_ 3 , _ 4 ^ 

Root [x3 — 3x2 _|_ 3, _ 4^ 3j ^ 

Root [19x3 — 19x2 — 3x — 1,2] —>• t2 
Root [38x3 _ 383.2 _|_ 403- _ 4^ 4j 
Root [38x3 _ 383.2 _|^ 403- _ 4^ 3j 
Root [x3 — 3x2 _|_ 3, _ 4^ 2] ^ A2 
Root [76x3 _ 703-2 _ 2a; — 1] —4^3 

Root [76x3 _ 703-2 _ 2x — 1,3] —>■ 713 
Root [x3 — x2 — 2x — 6,2] ^ 1/2 
Root [147x3 - 147x2 - 7x - 1,1] Ci 
Root [147x3 - 147x2 - 7x - 1,3] ^ Cs 
Root [588x3 _ 588a;2 + 77a; - 4,2] 02 

Root [1176x3 - 1176x2 - 161x - 6,1] pi 
Root [1176x3 - 1176x2 - 161x - 6,3] ^ pa 
Root [x3 — 3x2 _ 3. 4 2] —(T2 

Root [229x3 - 229x2 - 33x + 1,1] ^ n 
Root [229x3 - 229x2 - 33x + 1,3] T 3 
Root [229x3 - 229x2 + 5x + 2, 2] ^ 712 
Root [229x3 - 229x2 + 61x - 2,1]^ 

Root [229x3 - 229x2 + 61x - 2,3] C>3 

Root [x3 — 4x2 + 3x + 1,2] ^ X 2 
Root [49x3 — 49x2 + 1,1] —>• t/^i 
R oot [49x3 — 49x2 + 1,3] ^ 7/^3 
Root [49x3 — 49x2 — 14x + 1,2] -> a;2 
Root [49x3 — 49x2 + 14x — 1,1] —)■ Fi 
Root [49x3 — 49x2 + 14x — 1,3] —>■ Fa 
Root [x3 — 2x2 — 4x + 6,2] —>• 112 
Root [101x3 — 101x2 + 19x — 1,1] —>■ <7i 

Root [101x3 - 101x2 + 19x - 1,3] 73 

Root [202x3 - 202x2 - 42x - 1, 2] ^ £2 

Root [404x3 — 404x2 — 14x + 3,1] —>■ 

Root [404x3 — 404x2 — 14x + 3,3] —?> 7 I 3 

Root [x3 — x2 — 4x — 1, 2] xTa 
Root [169x3 — 169x2 — 26x + 1,1] —>■ wi 

Root [169x3 — 169x2 — 26x + 1,3] uj 3 

Root [169x3 — 169x2 — 13x + 5,2] —>• p 2 

Root [169x3 — 169x2 + 26x — 1,1] —>■ pi 

Root [169x3 — 169x2 + 26x — 1,3] —>■ tpa 


Table 6: Key to level sums table 
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11 Conclusion 


This paper has used a collection of multidimensional continued fractions to construct a fam¬ 
ily of sequences called TRIP-Stern sequences. These sequences reflect the properties of the 
multidimensional continued fractions from which they are generated. We have studied the 
sequences of maximal and minimal terms - and positions thereof. We have also character¬ 
ized the sums of levels and examined restrictions on terms appearing in a given TRIP-Stern 
sequence. We found that several of the level sum sequences or corresponding recurrence re¬ 
lations are well-known. Lastly, we introduced generalized TRIP-Stern sequences and proved 
several analogous results. 

We will conclude with a few unanswered questions: Do recurrence relations for row max¬ 
ima and their locations exist in general? We ask this because such relations could not always 
be found. What are the forbidden triples corresponding to (a, Tq, Ti) G Rf other than (e, e, e)? 
Do the distributions of terms in the TRIP-Stern sequences have any interesting properties? 
The terms of the TRIP-Stern sequences are the denominators of the convergents of the 
corresponding multidimensional continued fractions. Can these sequences reveal anything 
about approximation properties of multidimensional continued fractions? As discussed in 
Dasaratha et ah [10], many known multidimensional continued fractions are combinations of 
our family of 216 maps; as a result, it may be of interest to construct analogous sequences 
using select combination maps. 

There are polynomial analogs of Stern’s diatomic sequence (as in the work of Dilcher and 
Stolarsky [13, 14], of Coons [9], of Dilcher and Ericksen [12], of Klavzar, Milutinovic and 
Petr [30], of Ulas [50, 51], of Vargas [52], of Bundschuh [5], of Bundschun and Vaananen [6] 
and of Allouche and Mendes France [1]). What are the polynomial analogs for TRIP-Stern 
sequences? 

In essence, in this paper we start with a triple of numbers v = (a, b, c) and two 3x3 
matrices A and B and then examine the concatenation of the triples 

uA, vB, vAA, vAB, vBA, vBB, vAAA, vAAB, vABA ,.... 

Since there are 216 different triangle partition algorithms, we have 216 different pairs of 3 x 3 
matrices. Naively then, we would expect for there to be 216 different types of sequences, or, 
216 different stories. As we have seen, this is not the case. We have found clear patterns 
and classes among the 216 different TRIP-Stern sequences. The real question is why these 
patterns exist. Further, do the TRIP-Stern sequences that share, say common sequences of 
maximum terms, have common number theoretic properties? These questions strike us as 
hard. 
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